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ABSTRACT 
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usefulness  of  the  output  data  in  the  preliminary  design 
of  tracking  and  communications  networks. 


REVIEW  AI;D  approval 

This  technical  documentary  report  has  been  reviewed  and 
is  approved. 


VERNON  A.  DANDER 
Capt* ,  USAF 
Project  Officer 


iii 


TABLE  OF  CONTENTS 


Section  Pag;e 

1.  INTRODUCTION  1 

2.  USE  OF  VEHICLE  SIMULATION  IN  COMMAND  AND  CONTROL  STUDIES  1 

2.1  Line-of -Sight  Range  from  Each  Tracking  Site  to  the  2 

Re-entry  Vehicle 

2.2  Line-of -Sight  Range  Rate  2 

2.3  Elevation  Angle  of  Line-of-Sight  between  Each 

Tracking  Site  and  the  Vehicle  2 

2.4  Effect  of  Ground  Control  Delays  and  Errors  on 

Vehicle  Trajectory  3 

3.  COORDINATE  SYSTEMS  3 

3.1  Position  Coordinates  3 

3.2  Velocity  Coordinates  5 

3.3  Body  Attitude  Coordinates  5 

4.  ORGANIZATION  OF  COMPUTER  PROGRAMS  5 

5.  EQUATIONS  OF  MOTION  10 

6.  GUIDANCE  SYSTEM  10 

7.  ACCELERATIONS  AND  STRUCTURAL  TEMPERATURES  12 

8.  DOWNRANGE  AND  CROSSRANGE  PROJECTIONS  OF  VEHICLE  POSITION  13 

9.  TRACKING  STATION  MEASUREMENTS  14 

10.  RESULTS  OF  SIMULATIONS  16 

10.1  Analog  Results  16 

10.2  Digital  Results  25 

11.  CONCLUSIONS  30 


V 


APPEI®IX  A  --  DEAIVATIOU  OF  EQUATIOIJS  OF  MOTION  A-31 

APPEN'DIX  B  --  DEHIVATIOiJ  OF  CLOSED -FOHI-1  PREDICTION 

EQUATIONS  FOR  LIFTING  RE-ENTRY  TRAJECTORIES  -.-39 

appe:t)ix  c  --  guida::ce  systs:i  for  lifting  re-entry  'jeiilcu:  c-hT 

/J'PEI.'DIX  D  --  SYMBOLS  D-6l 

;>ppe:.dix  e  --  references  E-65 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 


4 

6 

7 

9 

15 

17 

18 

19 

20 

21 

22 

23 

24 

26 

27 


ILLUSTRATIONS 


Position  and  Velocity  Coordinate  Systems 
(1A-13,114) 

Downrange  and  Crossrange  Position  Coordinates 
(lA-13,113) 

Force  System  and  Vehicle  Body  Attitude 
(lA-13,116) 

Block  Diagram  of  Simulation 
(lA-13,117) 

Vehicle  Position  with  Respect  to  Tracking  Site 
(lA-13,115) 

Drag  Coefficient  Vs.  Lift-to-Drag  Ratio 
(lA-13,109) 

Angle  of  Attack  Vs.  Lift-to-Drag  Ratio 
(lA-13,108) 

Re-Entry  at  Constant  ~ 

(lB-13,105)  ° 

Typical  Near-Equilibrium  Glide  Trajectory 
(IB- 13, 106) 

Effects  of  Phugoid  Damping 
(lB-13,107) 

Re-Entry  Trajectories  for  Various  Bank  Angles 
(lA-13,315) 

Re-Entry  Trajectories  for  Various  Bank  Angles 
^  =  1.5  (lA-13,314) 

Re-Entry  Trajectories  for  Various  Bank  Angles 
^  =  1.0  (lA-13,313) 

Ground  Area  Attainable  at  Constant  Bank  Angle  a 

and  1.0  ^  ^  i  2.0  -  70°  i  a  ^  70°  (lA-13,312) 

Re-Entry  Vehicle  Trajectory 
(MA  2058) 


vii 


FIGURE 


PAGE 


16 

History  of  Altitude,  Velocity,  and  Stagnation  - 
Point  Skin  Temperature  (lA-13,327) 

28 

17 

Tracking  History 
(1A“13,328) 

29 

A-1 

Coordinate  Systems  for  Derivation  of  Equations 
of  Motion  (lA-13,311) 

A-32 

C-1 

Ground  Area  Attainable  (lA-13,111) 

C-l;9 

C-2 

Motion  of  Destination  of  Fixed  Reference  Sphere 
(lA-13,110) 

c-53 

C-3 

Relationship  Between  Commanded  “  Cos  c  and 

Normalized  Longitudinal  Error  (lA-13,112) 

C-58 

viii 


1 .  INTRODUCTION 

The  entry  of  ballistic  and  space  vehicles  into  the  earth *s  atmosphere 
has  long  been  of  interest  to  the  Air  Force.  Earlier,  interest  centered 
on  the  problems  associated  with  ballistic  missile  re-entry  vehicles. 

Today,  this  interest  has  been  broadened  to  include  manned  ballistic 
vehicles,  like  the  Gemini  and  Apollo  capsules,  and  maneuverable  lifting 
vehicles  of  the  X-20  type.  Accompanying  the  development  of  these  diverse 
vehicles,  there  has  of  course  been  an  impressive  growth  of  ground-based 
electronic  systems  supporting  vehicle  testing  and  operation.  As  the 
Air  Force’s  role  in  space  operations  develops,  a  continuing  need  can  be 
expected  for  quantitative  data  on  which  to  base  the  design  of  these 
electronic  systems. 

This  memorandum  describes  a  computer  simulation  of  re-entry  which 
has  been  developed  as  a  means  of  generating  some  of  the  quantitative 
information  needed  in  ground  system  design.  The  complete  simulation  is 
programmed  to  run  on  the  IBM  7030  digital  computer,  but  an  abridged 
version  has  also  been  developed  for  use  on  MITRE  * s  PACE  analog  computer. 

The  mathematical  model  which  forms  the  basis  of  the  computer  program  is 
discussed,  and  typical  results  produced  by  the  two  computers  are  presented. 

Section  2,  which  follows,  discusses  some  specific  questions  concern¬ 
ing  ground  system  design  for  which  the  simulation  can  provide  useful 
answers.  In  Sections  3  through  9,  the  organization  and  content  of  the 
simulation  are  described  and  important  assumptions  are  given.  Some 
actual  results  of  operating  the  simulated  system  on  both  analog  and 
digital  computers  are  presented  in  Section  10. 

2.  USE  OF  VEHICLE  SIMULATION  IN  COMMAND  AND  CONTROL  STUDIES 

The  actions  of  a  ground-based  electronic  system  supporting  the 
mission  of  an  aerospace  vehicle  are  intimately  related  to  the  motion  and 
activities  of  the  vehicle  itself.  Analysis  of  the  vehicle’s  mission  and 
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performance  is  therefore  a  prerequisite  for  determining  the  functional 
and  technical  requirements  of  such  a  ground  system.  Simulation  provides 
me  means  for  this  analysis  and  permits  interactions  between  the  vehicle 
and  the  ground  system  to  be  observed  while  these  individual  systems  are 
still  in  the  planning  stage^  Simulation  furnishes  the  quantitative 
information  necessary  for  specifying  quantitative  ground  system  require¬ 
ments.  In  the  case  of  the  present  simulation,  the  data  generated  and 
the  questions  these  data  can  help  to  answer  are  as  follows.  (All  data 
are  produced  as  functions  of  time.) 

2 . 1  Line-of -S ight  Range  from  Each  Tracking  Site  to  the 

Re-Entry  Vehicle 

A  time  history  of  range  between  the  vehicle  and  ground  stations 
is  a  necessity  in  the  rational  solution  of  range-dependent  tracking  and 
communication  problems.  Path  loss,  a  basic  variable  in  radar  tracking 
and  communication,  is  an  explicit  function  of  line-of-sight  range. 

Whether  the  immediate  design  problem  concerns  required  transmitter  power, 
communication  bandwidth,  or  the  contribution  of  receiver  thermal  noise 
to  radar  tracking  error,  range  is  a  variable  which  must  be  known. 

2 . 2  Line-of -Sight  Range  Rate 

Knowledge  of  the  range  rate  expected  at  a  given  ground  station 
during  the  tracking  period  is  useful  for  two  reasons.  First,  it  permits 
the  calculation  of  Doppler  shift  in  signals  sent  between  the  vehicle 
and  the  ground.  Secondly,  with  a  knowledge  of  a  given  radar's  range 
tracker  characteristics,  it  enables  one  to  estimate  the  range  tracking 
errors  that  would  be  produced  in  that  radar,  as  a  function  of  time. 

2 . 3  Elevation  Angle  of  Line  of  Sight  between  Each  Tracking 

Site  and  the  Vehicle 

By  observing  the  time  history  of  elevation  angle  at  each 
tracking  site,  it  is  a  simple  matter  to  determine  when  the  vehicle  is 
above  the  radar  horizon  at  that  station  and  thus  find  out  when  that 
station  can  track  and  communicate  with  the  vehicle,  given  that  the 
range  between  the  two  is  not  restrictive.  With  this  basic  information 
for  all  ground  stations,  and  with  a  knowledge  of  the  communications 
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traffic  through  each  station  when  it  is  actively  tracking,  a  time  history 
of  total  traffic  throughout  the  ground  communications  network  can  be 
constructed.  This  information  will  give  directly  the  traffic  handling 
requirements  of  the  various  segments  of  the  ground  network. 

The  history  of  vehicle  visibility  from  the  tracking  stations 
can  be  used  in  conjunction  with  the  corresponding  range  histories  to 
evaluate  the  effectiveness  of  various  deployments  of  ground  stations. 

In  comparing  various  deployments  of  a  fixed  number  of  stations,  the 
criterion  of  effectiveness  might  be  the  total  length  of  time  during  which 
the  vehicle  is  simultaneously  visible  (at  positive  elevation  angle)  and 
within  a  specified  range  of  one  or  more  stations.  Alternatively,  the 
problem  might  be  to  choose  a  re*'entry  trajectory  which  would  make  maxi-^ 
mum  use  of  an  existing  ground  system  configuration. 

2 . 4  Effect  of  Ground  Control  Delays  and  Errors  on  Vehicle 

Trajectory 

The  response  of  maneuverable  re-entry  vehicles  to  guidance 
and  control  actions  is  manifested  in  the  character  of  the  trajectory, 
in  the  time  histories  of  vehicle  accelerations  and  skin  temperatures, 
and  in  the  location  of  the  terminal  point  relative  to  the  intended 
destination.  The  effects  of  errors  and  time  delays  in  vehicle  guidance 
commands  can  be  determined,  therefore,  by  observing  these  same  quantities 
Ascertaining  these  effects  is  vital  in  any  study  aimed  at  determining 
the  feasibility  of  emergency  guidance  of  maneuverable  vehicles  by  the 
ground  system. 

3 .  COORDINATE  SYSTEMS 

It  will  be  useful  to  define  at  the  outset  the  coordinate  systems 
in  which  vehicle  position^,  velocity,  and  body  attitude  are  measured. 

3 . 1  Position  Coordinates 

Vehicle  position  is  measured  with  respect  to  a  rotating, 
spherical  earth.  The  three  position  coordinates,  shown  in  Fig.  1,  are 
longitude  0,  latitude  cp,  and  radius  r  from  the  center  of  the  earth. 
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A  secondary  set  of  coordinates,  the  downrange  and  cross-range 
distances  of  the  projected  position  on  the  earth*s  surface  from  the 
initial  position,  is  illustrated  in  Fig.  2.  The  downrange  axis  is  a 
great  circle,  the  plane  of  which  contains  the  vehicle *s  initial  velocity 
vector.  The  crossrange  axis  is  the  varying  great  circle  which  passes 
through  the  instantaneous  projected  position  and  intersects  the  downrange 
axis  at  right  angles.  Referring  to  Fig.  2,  downrange  and  crossrange 
distances  are  defined  as  the  length  of  the  arc  ^  along  the  downrange 
axis  and  the  length  of  the  arc  b  along  the  crossrange  axis,  respectively. 

3 . 2  Velocity  Coordinates 

Considered  as  a  point  mass,  the  vehicle  has  a  velocity  of 
magnitude  V  relative  to  earth  axes,  in  a  direction  defined  by  the 
angles  T)  and  y,  as  shown  in  Fig.  1.  The  angle  T)  represents  the  vehicle'? 
heading  with  respect  to  north;  the  angle  y  is  the  angle  of  the  velocity 
vector  above  the  local  horizontal  plane. 

3 . 3  Body  Attitude  Coordinates 

The  vehicle  is  assumed  to  maneuver  with  zero  sideslip  at  all 
times.  Body  attitude  is  then  defined  by  angle  of  attack  (y  and  bank 
angle  a,  as  depicted  in  Fig.  3.  Angle  of  attack,  of  course,  is  the 
angle  between  the  vehicle's  longitudinal  reference  axis  and  the  velocity 
vector.  The  bank  angle  is  the  angle  between  the  lift  vector  and  the 
vertical  plane  containing  the  velocity  vector. 

4.  ORGANIZATION  OF  COMPUTER  PROGRAMS 

A  re-entry  vehicle  may  be  any  of  a  number  of  types.  In  its 
simplest  form,  it  is  an  uncontrolled  body,  falling  ba llis tically  through 
the  atmosphere.  In  a  more  sophisticated  form,  it  may  contain  a  closed- 
loop  guidance  and  control  system  which,  by  modulating  the  aerodynamic 
forces  on  the  vehicle,  brings  it  to  some  specified  terminal  conditions 
of  position  and  velocity.  The  Mercury  capsule  is  an  example  of  the 
first  type,  whereas  the  latter  type  is  represented  by  the  X-20. 
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To  make  the  computer  programs  applicable  to  the  widest  range  of 
re-entry  vehicle  types  and  to  simplify  later  expansion  in  the  scope 
of  the  simulations,  the  programs  were  constructed  in  modular  form. 

In  the  case  of  the  digital  program,  this  means  that  nearly  all  compu¬ 
tations  are  done  in  subprograms,  and  the  main  program  is  in  effect  an 
executive  program.  Changes  in  the  system  being  simulated  are  easily 
made  by  changing  appropriate  subprograms. 

Figure  4  is  a  block  diagram  showing  the  over-all  organization  of 
the  simulation.  The  most  basic  part  is  that  identified  as  ’’Vehicle 
Dynamics”  —  the  subprogram  relating  position  and  velocity  of  the  vehicle 
to  bank  angle  and  aerodynamic  lift  and  drag  coefficients  Position 
and  velocity  are  derived  in  this  section  of  the  simulation  by  machine 
integration  of  the  differential  equations  of  motion.  Guidance  compu¬ 
tations,  based  primarily  on  instantaneous  position  and  velocity,  are 
performed  in  an  independent  subprogram.  The  guidance  system  determines 
the  desired  instantaneous  bank  angle  and  lift-to-drag  ratio.  In  the 
Aerodynamic  Data  subprogram,  lift-to-drag  ratio  is  translated  into 
angle  of  attack,  lift  coefficient,  and  drag  coefficient,  by  means  of 
tabulated  aerodynamic  characteristics  of  the  vehicle.  Local  atmospheric 
density,  required  in  many  computations,  is  computed  from  vehicle  altitude 
in  the  subprogram  called  ’’Model  Atmosphere”.  Operating  together,  these 
four  subprograms  produce  a  time  history  of  the  vehicie*s  position, 
velocity,  and  attitude. 

Three  other  subprograms  draw  upon  these  results  and  periodically 
compute  the  following  information: 

(a)  Downrange  and  crossrange  components  of  the  vehicle’s 
pos it ion . 

(b)  Specific  force  per  unit  mass  on  the  vehicle  (or 
’’acceleration”  felt  by  the  crew  in  a  manned  vehicle)  . 

(c)  Critical  structural  temperatures. 

(d)  Antenna  gimbal  angles  and  geometric  range  and  range  rate 
of  ground  based  stations  tracking  the  vehicle. 
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5.  EQUATIONS  OF  MOTION 

The  equations  of  motion  used  in  this  simulation  are  those  of  a 
point  mass  under  the  influence  of  lift,  drag,  and  gravitational  forces. 
Rotational  dynamics  of  the  vehicle  are  not  included. 

The  assumption  justifying  the  exclusion  of  rotational  dynamics  is 
that  if  the  guidance  system  calls  for  a  given  angle  of  attack  and  bank 
angle,  the  attitude  control  system  or  flight  control  system  of  the  vehicle 
can  produce  those  angles  with  a  negligible  time  lag.  Since  rotational 
time  constants  are  extremely  short  relative  to  the  total  time  of  flight, 
and  only  the  resultant  trajectory  of  the  vehicle  is  of  interest  here, 
this  assumption  introduces  no  significant  error  in  the  results. 

An  inverse'-square  gravitational  force  field  is  used  in  the  simula¬ 
tion.  To  first-order  accuracy,  therefore,  the  equations  ace  valid  for 
unpowered  vehicles  at  any  altitude,  including  satellites  out  of  the 
sensible  atmosphere,  vehicles  within  the  atmosphere,  and  re-entry  vehicles 
passing  from  free  space  through  the  atmosphere  to  the  surface  of  the 
earth . 

The  equations  of  motion  and  their  derivation  are  given  in 
Appendix  A. 

6.  GUIDANCE  SYSTEM 

Whereas  the  equations  of  motion  are  the  same  for  any  vehlcl-, 
there  is  virtually  no  limit  to  the  number  of  guidance  laws  which  may 
control  a  re-entry  vehicle's  trajectory.  For  the  case  of  lifting  re¬ 
entry  vehicles,  the  manifold  possibilities  in  guidance  system  design 
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are  enumerated  and  discussed  in  broad  perspective  in  Ref.  1.  As  the 
98-reference  bibliography  of  Ref.  1  attests,  there  is  a  large  body  of 
literature  on  the  subject  of  re-entry  guidance. 

The  computer  programs  are  so  organized  that  a  wide  variety  of 
guidance  techniques  can  be  simulated.  By  proper  choice  of  guidance 
equations  and  aerodynamic  data,  both  unguided  ballistic  re-entry 
bodies  and  maneuverable  re-entry  gliders  with  relatively  complex 
guidance  systems  can  be  accommodated. 

The  guidance  systems  used  in  producing  the  results  reported  in 
this  memorandum  are  discussed  in  Appendix  C.  The  analog  simulation 
used  a  constant  bank  angle  and  a  lift-to-drag  ratio  which  was  varied 
about  a  fixed  nominal  value  by  an  altitude  -  rate  damping  law.  The 
digital  simulation  contained  a  more  complex  scheme.  Guidance  commands 
were  based  on  both  altitude  rate  and  the  position  of  the  desired  destina¬ 
tion  with  respect  to  the  locus  of  attainable  destinations,  the  latter 
being  calculated  from  approximate,  closed-form  equations.  This  predic¬ 
tion  -  type  guidance  scheme  is  similar  to  those  described  in  Refs.  2, 

3,  and  4.  The  closed-form  prediction  equations  are  developed  in 
Appendix  B. 

It  has  often  been  suggested  that  a  ground-based  guidance  capability 
be  provided  for  maneuverable  re-entry  vehicles,  as  a  means  of  controlling 
the  trajectories  of  vehicles  which  have  suffered  primary  guidance  system 
failures.  Such  a  back-up  guidance  system  would  make  radar  measurements 
of  the  vehicle's  trajectory,  collect  telemetered  data  on  various  flight 
variables,  and  compute  guidance  commands  on  the  ground.  The  commands 
would  be  transmitted  to  the  vehicle  through  a  ground- to-vehic le  data 
link.  Many  problems  can  be  visualized  in  this  mode  of  guidance;  one  of 
these  is  the  possible  adverse  effect  of  data  processing  delays  on  the 
vehicle  trajectory.  To  permit  investigation  of  this  problem,  the  digital 
simulation  incorporates  the  option  of  introducing  a  time  delay  of  as 
much  as  60  seconds  between  calculation  of  the  guidance  commands  and 
application  of  these  commands  to  the  vehicle's  flight  control  system. 
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7  ACCELERATIONS  AND  STRUCTURAL  TEMPERATURES 

Accelerations  felt  by  human  crew  members  and  temperatures  of 
critical  parts  of  the  vehicle  skin  must  be  kept  under  certain  limits 
during  re-entry.  It  is  therefore  necessary  to  know  the  history  of 
these  variables  throughout  a  simulated  flight.  With  nomenclature  as 
defined  in  the  list  of  symbols,  components  of  specific  force  per  unit 
mass  along  the  vehicle's  body  axes  are  as  follows: 


Longitudinally  (positive  aft): 


^  [-D 


C_  Cos  CK  -  C  Sin  cy 


(1) 


Laterally,  in  the  plane  formed  by  the  longitudinal  axis 
and  the  lift  vector  (positive  upward): 


n  =  ^ Sin  (7  4*  C  Cos  (y 
n  mg  i_  D  L, 

e 


The  resultant  value  is 


n  = 


mg 


+  C 


2 

L 


(2) 


(3) 


In  this  simulation,  temperatures  are  computed  at  two  critical 
points  on  the  surface  of  a  winged  re-entry  vehicle:  the  stagnation 
point  and  the  bottom  of  the  wing.  The  equations  are  based  on  the 
assumptions  that  these  outer  surfaces  are  radiation-cooled  and  that 
the  flow  over  the  surfaces  is  laminar.  Stagnation-point  temperature 
is  given  by  the  equation  (Refs.  2  and  5). 


13880 

„  1/8 


/"lb.  sec.  \ 

V  ftX  /  -i 


1/8 


0.788 


2.3769  X  10 


]  R] 


degrees  F.  abs .  (4) 
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where 


(5) 


Wing-bottom  temperature  is  as  follows  (Ref.  2): 


"WB  ■  \b 


0.2  „0.467 


degrees  F .  abs . 


(6) 


8.  DOWNRANGE  AND  CROSSRANGE  PROJECTIONS  OF  VEHICLE  POSITION 

Downrange  and  crossrange  components  of  the  vehicle’s  projected 
position  have  been  defined  in  Section  3.1  and  in  Fig.  2.  Equations  for 
these  position  components  can  be  derived  from  Figs.  1  and  2  by  the 
standard  methods  of  spherical  trigonometry.  It  is  easily  verified  that 
calculation  of  these  distances  requires  the  solution  of  the  following 
set  of  equations: 


Cos  cp  Sin  (0"0Q- 


(8) 


Sin  c 


(9) 


(10) 


Downrange  distance  =  R^a 


(11) 


Crossrange  distance  =  R^b 


(12) 
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9.  TRACKING  STATION  MEASUREMENTS 

Ground-based  tracking  stations  are  represented  in  the  digital 

simulation.  For  each  of  ten  stations,  the  program  computes  the  range, 

range  rate,  azimuth  angle,  and  elevation  angle  that  could  be  measured 

bv  a  radar  tracking  the  re-entry  vehicle. 

Figure  5  shows  the  re-entry  vehicle’s  position  relative  to  the 

t  h 

coordinate  system  center  at  the  i  tracking  site.  In  terms  of  the 
local  coordinate  system,  vehicle  position  is 


X  . 
1 


=  r  Cos  cp  Sin  (0-0^) 


(13) 


y.  =  r  Sin  cp  Cos  cp^  ■  Cos  cp  Sin  cp^  Cos 


z  . 
1 


=  r  Sin  cp  Sin  cp.  ■+*  Cos  cp  Cos  cp.  Cos 


where 


r.  =  R  h. 
L  e  1 


=  R  -f  h 


(16) 


The  slant  range  from  the  tracking  station  to  the  vehicle  is 


(17) 


Azimuth  and  elevation  angles  are  given  by  the  equations 


(18) 


(19) 


L 


Ik 


lA-  13  115 


LOCAL 

VERTICAL 


VEHICLE 
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VEHICLE  POSITION  WITH  RESPECT  TO  TRACKING  SITE 

FIG,  5 
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By  differentiating  the  above  equation  for  and  using  Eqs .  13, 
lA  ,  and  15  above  and  A-21,  A-22,  and  A-23  of  Appendix  A  to  simplify  the 
result  ,  it  can  be  verified  that  range  rate  is  given  by  the  expression 


r  .  +  z  . 

1_ L 

r 


J 


-  Cosy 


Cos'll 


^uoscpb  incp^- 


(20) 


10.  RESULTS  OF  SIMULATIONS 
10 . 1  Analog  Results 

In  Figs.  8  through  13  there  are  shown  some  trajectories  pro¬ 
duced  by  the  analog  computer.  All  were  developed  for  a  lifting  re-entry 
vehicle  of  the  X-20  type,  with  the  aerodynamic  characteristics  given  in 
Figs.  6  and  7.  Figure  8  is  a  family  of  trajectories  at  constant  angle 
of  attack  (and  constant  lift-to-drag  ratio,  ~) ,  zero  bank  angle,  and 
varying  initial  flight  path  angle.  The  characteristic  lightly-damped 
phugoid  oscillation  is  apparent,  and  the  figure  shows  clearly  the  large 
effect  that  small  variations  in  initial  conditions  can  have  on  the 
motion  of  a  vehicle  of  this  type  if  re-entry  is  made  at  constant  L/D. 

In  Fig.  9,  altitude,  velocity,  dynamic  pressure,  stagnation- 
point  skin  temperature,  and  resultant  specific  force  per  unit  mass  during 
a  near-equilibrium  glide  are  presented  as  functions  of  time.  In  this 
case,  the  phugoid  motion  is  damped  by  an  altitude-rate  feedback  scheme 
which  varies  L/ D  around  a  nominal  value  of  1.5. 

The  practical  effects  of  damping  the  natural  trajectory 
oscillations  is  illustrated  in  Fig.  10.  Altitude,  s tagnat ion- point 
skin  temperature,  and  ’’acceleration"  are  presented  as  functions  of  time 
for  two  re-entry  trajectories  having  identical  initial  conditions.  In 
one  trajectory,  the  re-entry  vehicle  maintains  constant  L/D,  and  the 
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RE-ENTRY  TRAJECTORIES  FOR  VARIOUS  BANK  ANGLES 
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RE-ENTRY  TRAJECTORIES  FOR  VARIOUS  BANK  ANGLES 
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RE-ENTRY  TRAJECTORIES  FOR  VARIOUS  BANK  ANGLES 


natural  phugoid  oscillation  occurs;  in  the  second  trajectory,  the  phugoid 
damping  ratio  is  increased  to  around  0.7  by  altitude-rate  feedback.  The 
undesirable  temperature  and  acceleration  peaks  exhibited  along  the 
uncontrolled  trajectory  are  considerably  diminished  in  the  second  case. 

The  trajectories  presented  so  far  represent  flights  at  zero 
bank  angle.  All  lie  essentially  in  a  vertical  plane.  Figures  11,  12, 
and  13,  on  the  other  hand,  illustrate  the  lateral  maneuvering  capability 
of  a  lifting  re-entry  vehicle  when  non-zero  bank  angles  are  used.  In 
each  case,  these  figures  show  projections  of  vehicle  trajectories  on  the 
earth *s  surface  in  terms  of  latitude  and  longitude.  (Alternatively,  each 
line  is  the  locus  of  the  geocentric  sub-vehicle  point  during  re-entry.) 
Each  trajectory  was  generated  at  constant  bank  angle,  with  the  value 
indicated.  The  trajectories  represented  in  Fig.  11  were  made  with  a 
constant  L/ D  of  2.0.  Similarly,  Figs.  12  and  13  present  trajectories 
at  L/D’s  of  1.5  and  1.0,  respectively. 

The  vehicle  with  the  maneuvering  capability  depicted  in  these 
three  figures  can  reach  a  large  area  on  the  earth’s  surface.  If,  during 
re-entry,  bank  angle  is  restricted  to  values  between  -70^  and  -I-  70^,  and 
L/D  is  maintained  between  1.0  and  2.0,  the  "ground  area  attainable"  (GAA) 
existing  when  the  vehicle  is  in  its  initial  state  can  be  estimated  by 
superimposing  Figs.  11,  12,  and  13  and  drawing  the  envelope  of  the  points 
encompassed.  This  has  been  done,  and  the  result  is  shown  in  Fig.  14. 

By  the  proper  choice  of  bank  angle  and  L/D,  the  vehicle  is  able  to  reach 
any  point  inside  this  envelope. 

10 . 2  Digital  Results 

Figures  15,  16,  and  17  illustrate  the  kind  of  information  the 
digital  simulation  produces.  These  figures  present  data  on  a  single 
flight  of  a  lifting  re-entry  vehicle.  Aerodynamic  characteristics  of 
this  vehicle  were  identical  to  those  of  the  vehicle  for  which  the  analog 
data  were  derived,  but  in  this  case  the  vehicle’s  path  was  controlled  by 
the  re-entry  guidance  system  described  in  Appendix  C.  With  this  guidance 
system,  the  vehicle  was  able  to  reach  Edwards  AFB ,  California,  from  a 
starting  point  over  the  Indian  Ocean. 
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TRACKING  HISTORY 


The  horizontal  projection  of  the  vehicle*s  trajectory  on  the 
earth  is  shown  in  Fig.  15.  Time  markers  along  the  trajectory  permit 
correlation  of  the  vehicle*s  geographical  position  with  the  altitude, 
velocity,  and  skin  temperature  plotted  in  Fig.  16.  Also  shown  in  Fig. 

15  are  the  locations  of  eight  selected  tracking  stations  along  the 
vehicle*s  route.  In  this  simulated  flight,  it  was  postulated  that  each 
station  could  track  the  vehicle  while  the  latter  was  at  an  elevation 
angle  of  4  degrees  or  more.  The  resulting  tracking  coverage  throughout 
the  flight  is  illustrated  in  Fig.  17.  The  figure  presents,  for  each 
station,  a  shaded  bar  indicating  the  period  when  the  station  is  able  to 
track.  The  time  scale  is  the  same  one  used  in  Figs.  15  and  16.  Fig. 

17  shows  at  a  glance  when  the  tracking  network  as  a  whole  is  able  to 
collect  data,  and  how  many  stations  are  producing  data  simultaneously 
at  any  time. 

Figures  15,  16,  and  17,  taken  together,  should  give  the 
ground  system  designer  a  good  feel  for  the  fundamental  events  occurring 
during  this  flight,  and  provide  data  for  quantitative  analysis  of  system 
requirements.  Of  course,  information  on  vehicle  ”g^’  loads,  ground  area 
attainable,  tracking  station  range  rates,  and  other  variables  is  available 
also,  although  not  included  in  the  figures. 

11.  CONCLUSIONS 

This  memorandum  has  described  a  computer  simulation  of  re-entry 
which  has  been  developed  for  use  in  the  study  of  ground  electronic 
system  requirements.  Because  of  its  great  flexibility,  it  can  be  made 
to  represent  the  behavior  of  many  different  types  of  re-entry  vehicles. 

An  existing  MITRE  Corp.  simulation,  the  Flight  Trajectory  Profile 
program,  produces  trajectories  of  space  vehicles  throughout  the  boost 
and  exo- atmospheric  phases  of  flight.  Integrating  that  program  and 
the  present  re-entry  simulation  would  yield  a  useful  single  program 
capable  of  simulating  complete  space  missions,  from  initial  booster 
lift-off  to  final  landing  of  the  re-entry  vehicle. 
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APPENDIX  A 


Derivation  of  Equations  of  Motion 


BASIC  PRINCIPLE 

In  vector  form,  Newton’s  Second  Law  of  Motion  applied  to  a  point 
mass  m  is 

dv 

F  =  m  ,  (A-1) 

where  F  represents  resultant  force  and  V  is  velocity  measured  with 

Pl 

respect  to  inertial  space.  The  required  equations  of  motion  follow 
directly  from  this  relationship  when  proper  coordinate  systems  are 
defined . 

COORDINATE  SYSTEMS 

Coordinate  systems  have  already  been  defined  in  the  main  body  of 
this  report.  For  the  purpose  of  this  derivation,  however,  it  will  be 
convenient  to  use  two  additional  sets  of  orthogonal  reference  axes. 

These  axes,  an  shown  in  Fig.  A-1.  As 

the  figure  shows,  the  x^y^  plane  is  the  local  horizontal  plane;  x^  points 
northerly  and  z^  points  downward.  The  ^2^2^2  oriented  with 

respect  to  the  relative  velocity  vector  V;  X2  coincides  with  V  while 
y2  lies  in  the  local  horizontal  plane.  Both  sets  of  axes  are  right- 
handed  . 
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COORDINATE  SYSTEMS  FOR  DERIVATION  OF 
EQUATIONS  OF  MOTION 

FIG.  A-l 
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FORCES 


The  vehicle  moves  under  the  influence  of  three  forces:  gravity, 
aerodynamic  lift,  and  aerodynamic  drag.  The  latter  two  forces  are  defined 
with  respect  to  wind  axes;  that  is,  with  respect  to  axes  aligned  along 
and  perpendicular  to  the  relative  velocity  vector  V.  It  is  therefore 
convenient  to  expand  Eq.  A-1  in  terms  of  components  in  the 
sys  tern. 

The  drag  force  D  acts  opposite  to  V,  or  along  the  negative  x^  axis. 
Vehicle  banking  occurs  around  the  x^  axis,  and  the  lift  force  L  is  per¬ 
pendicular  to  V.  Therefore,  L  has  components  in  both  the  y^  and 
directions.  In  the  simple  central-force-field  model  of  the  earth  used 
here,  the  gravitational  force  is  directed  at  the  center  of  the  earth 
coordinate  system  and  has  components  in  the  x^  and  directions.  The 
resultant  force  vector  is  the  following: 


2y2"2 


D  -  mg  Sin  Y 
L  Sin  a 

L  Cos  a  ■+■  mg  Cos  Y 


(A- 2) 


The  subscript  that  the  components  given  are  measured  with 

respect  to  the  X2y222  Qxis  system. 

Lift  and  drag  forces  depend  on  dynamic  pressure  q  and  the  corres¬ 
ponding  aerodynamic  coefficients: 


L  =  q  A  C 

L 

(A-3) 

D  =  q  A 

(A-4) 
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where 


q  =  2  P  ^ 


(A-5) 


and  the  coefficients  C  and  C  are  taken  to  be  functions  of  angle  of 

Li  JJ 

attack  only.  At  the  hypersonic  velocities  existing  throughout  most  of 
the  re-entry  period,  neglecting  the  dependence  of  C  and  C  on  Mach 

Li  D 

number  does  not  introduce  serious  error. 

Gravitational  acceleration  g  is  determined  by  the  inverse-square 

law , 


g 


GM 

2 


(A-6) 


where  GM  is  the  gravitational  parameter  of  the  earth. 


ACCELERATIONS 

Like  the  forces,  the  rate  of  change  of  V  on  the  right-hand  side 

A 

of  Eq.  A-1  will  be  expanded  into  components  in  the  system. 

This  axis  system  moves  with  the  vehicle,  and  has  both  translational  and 
rotational  motion.  Because  of  the  angular  velocity  of  the  reference 
axes,  the  total  rate  of  change  of  V  is  as  follows: 

A 


TOTAL 


A- 


Ldt 


+  uo  X  V, 


(A- 7) 


The  first  term  on  the  right-hand  side  is  the  rate  of  change  of  V  rela- 

A 

tive  to  the  rotating  coordinates.  The  second  term  is  the  cross  product 
of  the  absolute  angular  velocity  of  the  coordinate  system,  with  V  . 

A 

Both  terms  are  vectors  and,  as  already  stated,  are  to  be  evaluated  with 
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respect  to  the  system.  Reference  6,  among  others, 

gives  the  theoretical  basis  for  Eq.  A-7. 

The  absolute  velocity  can  be  expressed  in  terms  of  the  relative 
velocity  V  and  earth *s  rate  Q  as  follows: 


"ax' 

r  • 

V  +  fjr  Cos  cp  Sin  T]  Cos  y 

<C 

> 

= 

Qr  Cos  cp  Cos  T] 

Qr  Cos  cp  Sin  T]  Sin  y 

■  ^ 

(A-8) 


The  subscript  denotes  that  the  components  given  are  measured  in 

the  system.  The  right-hand  side  of  Eq.  A-8  is  a  column  matrix 

of  the  components  of  the  vector  . 

Angular  velocity  of  the  system,  measured  with 

respect  to  the  axes,  is 


x^yi^l 


P  Cos  cp  -  Y  Sin  T] 

-  cp  +  Y  Cos  T] 

-  p  Sin  cp  +  T] 


(A-9) 


It  follows  from  Fig.  A-1  that  the  corresponding  expression  for  this 
vector  in  ^2^2^2  coordinates  can  be  found  by  applying  the  following 
trans  formation; 


^1^2  2 
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-  Sin  T] 
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Sin  T]  Cosy 
Cos  T] 

Sin  T|  SinY 


“Sin  Y 

0 

Cos  Y 


[=:] 


(A- 10) 
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This  transformation  leads  to  an  evaluation  of  fju  in  the  desired 


reference  frame: 
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The  cross  product  in  Eq.  A- 7  is  equivalent  to  the  following: 
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The  absolute  rate  of  change  of  V  can  now  be  found  by  using  Eqs . 

A 

A-8,  A“ll,  and  A-12  to  carry  out  the  operations  specified  in  Eq .  A- 7 , 
The  lengthly  results  thus  obtained  can  be  simplified  considerably  by 
the  use  of  the  following  identities: 


r  0  Cos  cp  Sin  Tj  +  r  9  Cos  T]  =  V  Cos  y 


(A-13) 


r  6  Cos  9  Cos  T]  -  r  9  Sin  T]  =  0 

r  =  V  Sin  y 


(A- 14) 
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The  angular  rate  0  can  be  eliminated  from  the  equations  by  noting 
that  it  is  equal  to  0  plus  the  earth's  inertial  rotational  rate: 


The  desired  result  is 


V2"2 


V  +  0  r  Cos  cp  [sin  cp  Cos  T]  Cos  y  -  Cos  cp  Sin  yj 


VCosy[T|-  Tan  cp  Sin  T|  Cos  y] 

+  2r}V[Cos  cp  Cos  T)  Sin  y  -  Sin  cp  Cos  y] 
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FINAL  EQUATIONS  OF  MOTION 

Three  final  equations  of  motion  result  when  Eqs .  A-2  and  A-17  are 
substituted  into  Eq.  A-1  and  corresponding  components  of  this  vector 
equation  are  set  equal  to  each  other: 
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V  =  -  “  -  g  Sin  y  -  n  r  Cos  p  [SinpCosTjCosy-CospSiny  ]  (A-18) 
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Since  these  are  first-order  differential  equations,  only  three  integrations 
are  involved.  Three  more  integrations  are  required  before  the  complete 
solution  can  be  obtained.  These  are  provided  by  the  following  kinematic 
equat ions : 


r  =  V  Sin  y 
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r  Cos  cp 
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Eqs.  A-18  through  A-23  constitute  the  complete  set  of  equations  of  motion. 

Lastly,  it  is  evident  from  the  geometry  of  the  problem  that  vehicle 
altitude  is  determined  by  the  relation 


h  =  r  -  R 
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APPENDIX  B 


Deviation  of  Closed-Form  Prediction  Equations  for 

Lifting  Re’-entry  Trajectories 

EQUILIBRIUM  GLIDE  TRAJECTORY 

The  closed-form  trajectory  prediction  equations  to  be  developed 
here  are  based  on  the  so-called  ’’equilibrium  glide”.  This  type  of 
trajectory  is  defined  as  one  along  which  the  centrifugal  and  vertical 
lift  forces  on  the  vehicle  exactly  counterbalance  the  gravitational 
force.  One  of  the  principal  merits  of  the  equilibrium  glide  is  that 
it  can  be  described  by  relatively  simple  analytical  relationships. 


HORIZONTAL  RANGE  PREDICTION 

If  two  of  the  equations  of  motion  derived  in  Appendix  A  (Eqs.  A- 18 
and  A-20)  are  written  in  terms  of  coordinates  fixed  to  a  non-rotating 
earth  model,  they  become  the  following: 


V.  =  -  -  -  g  Sin  Y. 
A  m  °  ^A 
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where  the  relative  velocity  V  has  been  transformed  to  absolute  velocity 

V  ,  and  a  flight  path  angle  Ya  defined  as  the  angle  of  the  vector 

A  A 

V  above  the  local  horizontal.  The  bank  angle  a*  is  defined  with  res- 

A  A 

pect  to  a  vertical  plane  containing  V  .  Now,  the  equilibrium  glide  is 

Pi, 

the  path  for  which  -y  =  0.  That  is, 

Pi 


L  Cos  g  Cos  Cos  y^ 
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and  there  is  equilibrium  between  lift,  gravity,  and  centrifugal  forces.  If 

flight  path  angle  y  is  restricted  to  small  values,  so  that  Cos  y  s-  L , 

A  A 

this  equation  is  equivalent  to 


L  Cos  a, 
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The  quantity  gr  is  equal  to  ,  the 
satellite  in  circular  orbit  at  radius  r. 
vritten  as 


square  of  the  velocity  of  a 
Therefore,  Eq.  B‘4  can  be 
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Clearly,  the  ratio  —  must  be 
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L  cos  a  ,  is  to  act  in  the  upw 

Pi 

With  the  assumption  that 

g  Sin  V,  «  Eq.  B-1  can  be 
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less  than  1  if  the  vertical  lift  componen 
ard  direction. 

the  flight  conditions  are  such  that 
simplified  to  the  following: 
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This  equation  can  be  combined  with  B-5  to  produce  a  useful  range 
prediction  result  when  velocity  is  related  to  horizontal  distance 
traveled.  For  small  values  of  y  , 

Pi 


dt 


\  cos 
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Where  dS  is  the  differential  horizontal  distance  traveled  in  time  dt. 
It  follows  that 
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A  dt  dt  dS  A  dt 


When  combined,  equations  B~5,  B-6 ,  and  B-8  yield  the  result 
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This  different lal  equation  in  V  can  be  integrated  in  closed  form  if 
L  ^ 

“  Cos  a*,  g*  and  V  are  held  constant.  The  latter  two  variables  actually 
D  A  c 

vary  very  little  during  lifting  re-entry  of  the  earth’s  atmosphere.  The 
quantity  “  Cos  is  controllable,  and  can  be  held  fixed  arbitrarily. 

The  result  of  integrating  Eq.  B-9  between  the  velocity  limits  V  and 

A 

zero  and  the  distance  limits  zero  and  S  is 


(B-IO) 


where  r  =  —  =  Average  value  of  radius  during  flight.  The  minus 
AV .  g 

sign  is  consistent  with  the  requirement  that  initial  re-entry  velocity 
be  sub-circular,  so  that  0  <  —  <  1  and  the  logarithm  on  the  right  hand 
side  is  always  negative.  This  equation  gives  the  approximate  horizontal 
distance  traveled  by  a  lifting  re-entry  vehicle  on  an  equilibrium  glide 
trajectory  with  constant  ^  Cos  a. ,  from  the  point  at  which  inertial 
velocity  is  V  . 

In  the  special  case  for  which  a .  =  0 ,  approximate  horizontal  dis- 
tance  is 
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AV. 


2 

2  J 


(B-11) 


The  projection  of  this  trajectory  on  the  surface  of  the  earth  is  a  great 
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circle  which  subtends  an  angle  \  approximately  equal  to 
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Therefore , 
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ALTITUDE  RATE 

In  this  section,  a  simple  expression  for  altitude  rate  during  an 
equilibrium  glide  will  be  developed  from  the  basic  equations  of  motion 
and  an  analytical  model  of  the  atmosphere. 

Altitude  rate  can  best  be  evaluated  indirectly,  by  writing  it  as 
follows : 


dh 

dt 
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The  problem  thus  becomes  one  of  evaluating  the  derivations 
■  dV. 


dt 


and 


An  approximation  for  has  already  been  given  in  Eq.  B-6: 


dt  m 


(B-14) 


The  derivative  can  be  derived  from  the  basic  equation  of  the 

flight  path,  Eq.  B-5. 

The  aerodynamic  lift  force  L  in  this  equation  is  given  by  the 
expression 
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where  V  is  the  velocity  of  the  vehicle  relative  to  the  earth's  rotating 
atmosphere.  The  velocities  V  and  differ  by  a  maximum  of  about  1500 
feet  per  second  at  the  equator  and  less  at  other  latitudes ~a  relatively 
small  percentage  difference  during  most  of  the  re-entry  glide.  As  an 
approximation,  therefore,  can  be  used  for  V,  and  Eq.  B-5  becomes 
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By  taking  logarithms  of  both  sides  of  this  equation  and 

with  respect  to  altitude  h  while  holding  ^  Cos 

mg 


the  following  result  is  obtained: 


differentiating 

and  V  constant, 
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By  the  use  of  Eq.  B-5,  this  can  be  simplified  as  follows: 
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In  evaluating  the  right-hand  derivative  in  this  equation,  it  is  convenient 
to  use  the  exponential  model  atmosphere: 
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It  can  be  verified  that  by  virtue  of  this  relationship  between  density 
and  altitude,  the  following  is  true: 
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dV 

Therefore,  the  derivative  (Eq.  B-18)  becomes 
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With  the  substitution  of  this  derivative  and  the  one  of  Eq.  B- 14 
into  Eq.  B-13,  altitude  rate  is  found  to  be 
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or 
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PREDICTION  OF  REMAINING  TIME  OF  FLIGHT 

Equation  B-7  provides  a  simple  means  for  determining  the  time  of 
flight  remaining  at  any  instant: 

dt  (B-24) 
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By  combining  Eq.  B-9  with  this  one,  the  following  equation  in  t  and  V 

Pi 

can  be  obtained: 
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If  the  indicated 

V  held  constant 
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velocity  is  V. 

A 


integrations  are  carried  out,  with 
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Cos  cr. 
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,  the  time  of  flight  remaining  at  the  instant  when 
is  found  to  be 


(B-26) 
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APPENDIX  C 


Guidance  System  for  Lifting  Re-entry  Vehicle 


SCOPE 

This  appendix  describes  a  system  for  guiding  a  lifting  re-entry 
vehicle,  such  as  the  X-20  glider,  to  a  desired  destination  on  the  earth 
with  high  probability  of  success.  This  system  is  only  one  of  many  possi¬ 
bilities.  It  differs  from  the  one  actually  used  in  the  X-20,  but  is  very 
similar,  in  concept  and  in  many  details,  to  the  energy  management  system 
described  in  Ref.  2.  Basically,  it  is  designed  to  control  the  locus  of 
possible  landing  points. 

Re-entry  of  a  maneuverable  vehicle  like  the  X-20  is  a  process  invol¬ 
ving  several  distinct  phases.  While  the  vehicle  is  still  out  of  the 
atmosphere,  its  path  is  ballistic,  and  the  only  control  exercised  is 
that  of  vehicle  attitude.  When  aerodynamic  forces  develop,  the  trajec¬ 
tory  becomes  controllable  and  a  second  phase,  the  glide  phase,  begins. 

It  IS  during  this  phase  that  the  critical  problems  of  energy  management 
occur.  When  the  vehicle  nears  its  destination  and  most  of  its  initial 
energy  has  been  dissipated,  control  of  the  vehicle  passes  to  a  terminal 
guidance  system,  and  the  last  phase  begins.  This  transition  occurs  at 
an  altitude  of  approximately  100,000  feet,  where  velocity  is  on  the  order 
of  4000  feet  per  second.  The  terminal  guidance  system,  which  may  be  a 
human  pilot  using  visual  cues,  flight  instruments  and  advisory  information 
radioed  from  the  ground,  has  the  job  of  landing  the  vehicle  safely  in  a 


designated  recovery  area 


The  guidance  system  described  here  is  designed  to  carry  the  vehicle 
through  the  second  phase  of  re-entry.  Simulations  using  this  guidance 
•-vstem  ordinarily  end  when  an  altitude  of  100,000  feet  is  reached,  since 
this  mode,  of  flight  is  not  valid  beyond  that  point. 

GUIDANCE  CONCEPT 

Two  broad  categories  of  re-entry  guidance  systems  have  been  proposed 
and  studied  in  recent  years.  These  are  the  predic t ion- type  systems  and 
the  systems  designed  to  guide  the  vehicle  onto  a  pre-computed ,  reference 
trajectory-  The  present  scheme  belongs  in  the  first  category.  It  is 
based  on  the  concept  of  the  "ground  area  attainable"  (GAA)  or  "footprint" 
of  a  maneuverable  re-entry  vehicle.  The  GAA  is  defined  as  the  ground 
area  to  which  the  vehicle  is  capable  of  maneuvering  from  its  existing 
state  of  position  and  velocity.  In  Fig.  C-1,  the  GAA  is  shown  with  res¬ 
pect  to  the  re-entry  vehicle  in  a  plan  view  which  neglects  the  eartli's 
curvature.  The  intended  destination  on  the  ground  is  shown  at  an  arbitrary 
position  inside  the  GAA  boundary. 

As  the  vehicle  descends  through  the  atmosphere  and  its  velocity 
decreases,  the  GAA  shrinks  in  size  and  changes  shape  slightly.  In  order 
for  the  vehicle  to  reach  a  given  destination,  that  destination  must  be 
kept  within  the  diminishing  GAA.  From  an  alternate  point  of  view,  the 
GAA  must  be  controlled  by  vehicle  maneuvering  so  that  it  always  covers 
the  desired  destination. 

The  present  guidance  system  is  designed  to  drive  the  GAA  continuously 
to  the  position  in  which  the  desired  destination  is  exactly  in  the  center. 
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GROUND  AREA  ATTAINABLE 

FIG.  C-l 
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With  the  target  point  ideally  centered  in  the  GAA ,  random  perturbations 
of  the  trajectory,  produced  by  such  real  factors  as  guidance  equipment 
errors  and  vehicle  system  malfunctions  causing  temporary  loss  of  guidance 
or  control,  have  the  least  effect  on  the  ability  of  the  vehicle  to  reach 
that  point. 

Another  feature  of  this  guidance  scheme  is  its  use  of  the  equilibrium 
glide  trajectory  as  a  basis  for  computing  the  maneuvering  capability  of 
the  vehicle  at  any  instant.  The  GAA  at  any  time  is  determined  by  calcu¬ 
lating  the  limits  to  which  the  vehicle  could  travel  ^  equilibrium  glide 
tra lector les .  These  limits  are  not  the  absolute  limits  achievable  by 
optimum  programming  of  the  guidance  commands,  but  are  the  limits  which 
can  be  realized  with  the  present  guidance  system. 

The  two  basic  guidance  commands  are  lift-to-drag  ratio  (^)  and  bank 
angle  C  Both  are  based  on  the  computed  position  of  the  target  point  with 
respect  to  the  four  cardinal  points  of  the  GAA.  The  latter  points  are 
determined  by  means  of  approximate,  closed-form  equations.  Figure  C-1 
shows  the  cardinal  point  locations. 

GAA  CARDINAL  POINTS 

The  desired  landing  point  and  the  GAA  cardinal  points  are  specified 
with  respect  to  a  particular  reference  line  on  a  non-rotating  spherical 
earth.  This  line  is  the  great  circle  contained  in  the  plane  passed 
through  the  center  of  the  earth  and  the  vehicle's  instantaneous  velocity 
vector.  The  closest  and  most  distant  cardinal  points,  lying  on  the  refer¬ 
ence  great  circle,  are  reached  with  zero  bank  angle.  By  virtue  of  the 
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results  derived  in  Appendix  B  (Eq.  B-12) ,  equilibrium  glide  trajectories 
from  the  present  position  to  these  points  cover  the  following  central 
angles : 

Minimum  range: 


^MIN. 

Maximum  range: 

^MAX. 


(C-1) 
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where 
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and  V  is  the  velocity  of  the  vehicle  with  respect  to  inertial  space. 

In  terms  of  the  variables  used  in  the  equations  of  motion,  the  magnitude 

of  V.  is 
A 


-I 


V  +  2r2  rV  Cos  cp  Sin  T]  Cos  y  ■*"  (^r  Cos  cp)  (C-3) 


Similarly,  the  central  angle  from  each  lateral  cardinal  point  to  the 
reference  great  circle  is  approximately  as  follows: 


V  =  0.160 


1.19 

MAX. 
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This  last  relationship  was  derived  empirically  from  re-entry  trajectories 

produced  by  an  analog  computer.  It  is  valid  for  values  of  —  between 

1.0  and  2.0  and  for  values  of  —  between  0.2  and  0.9. 

c 

These  prediction  equations  (C-1,  C-2,  and  C-4)  are  reasonably 
accurate  when  applied  to  the  non-rotating  reference  frame  in  which  they 
were  derived. 


POSITION  OF  DESIRED  DESTINATION 

As  illustrated  in  Fig.  C-1,  the  position  of  the  destination  is 
defined  by  two  central  angles:  A>  the  one  between  the  destination  and 
the  vehicle,  and  |j, ,  the  one  subtended  by  the  great  circle  arc  passing 
through  the  destination  and  intersecting  the  reference  great  circle  at 
right  angles.  These  central  angles  are  defined  with  respect  to  a  non¬ 
rotating  reference  sphere  with  center  at  the  center  of  the  earth.  In 
Fig.  C-1,  the  position  of  the  destination  with  respect  to  this  reference 
system  is  the  position  it  occupies  at  the  end  of  the  vehicle*s  flight. 

By  analyzing  the  spherical  geometry  involved,  with  the  aid  of  Figs. 
C-1  and  C-2,  it  can  be  shown  that 


A  =  Cos  Sin  cp  Sin  cp  +  Cos  cp^  Cos  cp  Cos  (0 
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MOTION  OF  DESTINATION  ON  FIXED  REFERENCE  SPHERE 
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and 
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Tlie  angle  Ot  in  the  last  equation  is  the  angle  through  which  the 
r 

earth  rotates  in  the  remaining  time  of  flight,  t  .  The  equation  for 

F 

estimating  t  ,  derived  in  Appendix  B  (Eq.  B-26) ,  contains  the  factor 
L  ^  . 

—  cos  a  .  It  is  shown  in  a  later  section  of  this  appendix  that  the 
^  ^  r  L 

guidance  system  attempts  to  drive  this  quantity  to  the  value  h 

^D^MIN  J’  Therefore,  with  the  latter  expression  used  in  place  of  the 


■r  cos  c,  ,  Eq.  B-26  can  be  approximated  as 
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NORMALIZED  ERROR  SIGNALS 

Guidance  commands  are  computed  from  error  quantities  representing 
the  amount  by  which  the  destination  is  displaced  from  the  GAA  center. 
However,  to  match  the  magnitudes  of  the  commands  to  the  maneuvering 
capabilities  of  the  vehicle  at  all  times,  the  computed  errors  are 
normalized  When  the  destination  is  centered  in  the  GAA,  errors  are 
zero,  and  when  the  destination  coincides  with  one  of  the  cardinal  points 
the  magnitude  of  the  corresponding  error  is  unity.  Normalized  lateral 
error,  used  in  determining  the  required  bank  angle,  is  as  follows: 


Vc 
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Normalized  longitudinal  error,  used  to  compute  required  — ,  is  given  by 
the  following  equation: 


c-54 
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PHUGOID  DAMPING 


Bank  angle  command  is  based  on  normalized  errors  alone.  Lif o-drag 
ratio  command,  on  the  other  hand,  has  two  components;  one  based  on 


and  the  other  based  on  the  vehicle's  instantaneous  altitude  rate.  The 
latter  component  has  the  function  of  damping  out  the  phugoid  oscillation 
in  the  trajectory. 

The  phugoid  oscillation  is  a  lightly-damped  oscillatory  deviation 
from  a  smooth,  monotonic  trajectory.  The  nature  of  the  phenomenon  is 
illustrated  in  Fig.  8  of  the  main  body  of  this  paper.  If  the  flight 
conditions  of  a  re-entry  glider  flying  at  constant  ^  Cos  a  deviate  from 
those  for  an  equilibrium  glide  at  that  value  of  —  Cos  a,  a  phugoid 
oscillation  is  induced.  The  mechanism  of  the  oscillation  hinges  on  the 
variation  of  atmospheric  density  with  altitude,  and  is  discussed  in 
some  detail  in  Refs.  2  and  7. 

If  phugoid  oscillations  are  not  controlled,  vehicle  acceleration 
and  skin  temperatures  will  have  an  oscillatory  history  and  may  reach 
intolerably  high  peak  values.  It  is  therefore  desirable  to  damp  the 
oscillations  by  varying  ^  Cos  a,  and  avoid  high  peak  accelerations  and 
temperatures.  In  the  present  system,  the  required  damping  signal  is 
derived  from  the  difference  between  the  vehicle's  altitude  rate  and  that 
which  the  vehicle  would  have  if  it  were  on  a  reference  equilibrium  glide 
trajectory  with  its  present  velocity.  It  is  shown  in  Appendix  B  that 
this  equilibrium  altitude  rate  is 
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The  damping  signal  takes  the  form  of  an  increment  in  the  commanded  value 
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a  (  i  Cos  a  )  cr  (  ti  -  ) 


Reference  2  shows  that  a  more  nearly  constant  damping  ratio  can  be 
achieved  throughout  the  flight  by  making  the  proportionality  constant 
in  turn  proportional  to  vehicle  velocity.  Thus, 

A  (ices  a)-  (h  -  hj  (c-14) 


GUIDANCE  COMMANDS 

Bank  Angle 

Commanded  bank  angle  a  is  made  directly  proportional  to  the  normalized 
lateral  error,  with  a  proportionality  constant  weighted  by  the  normalized 
longitudinal  error: 
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This  angle  is  not  allowed  to  exceed  certain  limits,  however: 


^  a, 


MAX. 


(C-16) 


If  the  angle  computed  from  Eq.  C-15  falls  outside  this  allowable  range, 
the  value  of  the  nearest  limit  is  used  instead.  The  angle  determined 
by  Eqs.  C-15  and  C-16  is  used  in  the  subsequent  calculation  of  the  ~ 
command . 

Lift-to-Drag  Ratio 

As  Appendix  B  shows ,  longitudinal  range  is  determined  by  the  quan¬ 
tity  ^  Cos  a-  Lift-to-drag  ratio  command  is  computed  by  determining 
the  required  value  of  ^  Cos  a  and  dividing  this  by  Cos  a. 
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When 


The  first  of  the  two  components  of  —  Cos  a  is  based  on  ARp^. 

Vf  L  ^  ^ 

ARj^  is  zero,  ”  Cos  a  must  have  a  value  corresponding  to  the  landing 

point  halfway  between  the  two  longitudinal  cardinal  points: 
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MAX. 
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Non-zero  values  of  AR_  call  for  variations  in  —  Cos  a  from  this  central 

^  L  ^  ^ 

value.  The  function  used  to  relate  “  Cos  a  to  ARj^  is  given  in  Fig.  C-3. 

The  non-linearity  of  this  function  provides  an  over-control  feature. 

which  tends  to  center  the  GAA  over  the  destination  more  quickly  than 

* 

would  a  function  linear  throughout  the  range  -1  ^  ARp^  ^  +  1. 

This  component  of  “  Cos  a  is  called  the  nominal  value,  and  can  be 
expressed  analytically  as  follows: 

For  AR;  ^  -  AR;^: 

ft  “4om.  ■  ft4iK."°= " 
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For  0  <  AR*  <  AR*j^ 
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RELATIONSHIP  BETWEEN  COMMANDED  COSa 
AND  NORMALIZED  LONGITUDINAL  ERROR 

FIG.C-3 
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For  ARj^ 
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The  second  component  of  the  “  Cos  a  command  is  the  phugoid  damping 

term,  Eq.  C-14.  The  equilibrium  altitude  rate  h  used  in  computing 

e  L 

this  quantity  is  based  on  the  nominal  value  of  —  Cos  a  immediately 
above.  When  the  damping  term  is  expanded  by  the  use  of  Eq.  C-13  and  added 
to  the  nominal  value,  the  final  “  Cos  a  command  is  produced: 
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NOM. 


Lift-to-drag  ratio  is  developed  by  simply  dividing  the  last  equation  by 
Cos  Q. 

Finally,  the  command  generated  in  this  way  is  subjected  to  a  limiting 
process,  to  keep  the  lift-to-drag  ratio  of  the  vehicle  within  specified 
limits  at  all  times: 

"L  I  <  L  <  L 

If  the  value  of  “  computed  from  Eq. 
range  (Eq.  C-20)  ,  the  nearest  limiting 


C-19  falls  outside  the  allowable 
value  is  used  instead. 
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.APPEUDIX  D 

SYMBOLS 


Downrange  central  angle. 

Reference  area. 

Azimuth  angle  from  tracking  station  i  to  vehicle. 
Crossrange  central  angle. 

Central  angle  between  initial  and  present  positions. 
Drag  coefficient. 

Lift  coefficient. 

Aerodynamic  drag  force. 

Elevation  angle  of  line  of  sight  from  tracking 
station  i  to  vehicle. 

Resultant  force  on  vehicle. 

Vehicle  stagnation  point  design  parameter  (ft.) 

Local  gravitational  acceleration. 

Gravitational  acceleration  at  earth*s  surface. 
Universal  gravitationa 1  cons tant 
Vehicle  altitude  (positive  upward). 

Elevation  of  tracking  station  i  above  reference 
spherical  earth. 

Exponential  decay  parameter  of  earth’s  atmospheric 
dens ity . 

2  2 

Phugoid  damping  gain  constant  (sec.  /ft.  ). 


.Or 

Wing  bottom  temperature  constant  ^ 

(Ib.)^'^(ft.)^-^^^ 
Bank  angle  gain  constant  (degrees) . 

Bank  angle  cross-coupling  gain  constant  (degrees). 

Lift  force. 

Mass  of  vehicle. 

Mass  of  the  earth 

Resultant  specific  force  per  unit  mass  on  the 
vehicle  (in  units  of  earth ’ s -surf ace  g*s) . 
Longitudinal  component  of  specific  force  per  unit 
mass  (in  earth ’ s -surf ace  g’s). 


n 


N 


q 

r 
r . 

L 

R 

e 

R. 

1 


AR* 


V 

c 


X.  , 
L 


a 

3 


Y 


=  Lateral  component  of  specific  force  per  unit  mass 
(in  earth  * s -surface  g*s). 

=  Dynamic  pressure 

=  Vehicle*s  radius  from  earth*s  center 

=  Radius  from  earth*s  center  to  tracking  station  i. 

=  Radius  of  the  earth. 

=  Line-of -s ight  range  between  vehicle  and  tracking 
station  i 

=  Radius  of  vehicle  structure  at  stagnation  point  (ft.) 

=  Normalized  lateral  guidance  error 

=  Normalized  longitudinal  guidance  error 

=  Lift-to-drag  ratio  gain  parameter. 

=  Horizontal  distance  traveled  by  vehicle. 

=  Time 

=  Length  of  flight  time  remaining  at  any  instant. 

=  Stagnation-point  skin  temperature  (°F.  abs . ) 

=  Wing-bottom  temperature.  (^F.  abs.). 

=  Vehicle  velocity  with  respect  to  the  rotating  earth. 

=  Vehicle  velocity  with  respect  to  an  inertial 
reference  frame. 

=  Inertial  velocity  of  a  satellite  in  circular  orbit 
at  the  vehicle's  present  altitude  above  the  earth 
(=  gr)  . 

=  Vehicle  position  components  with  respect  to 

coordinate  system  centered  at  tracking  station  i. 

=  Angle  of  attack  (positive  upward) 

=  Angle  between  vehicle's  local  meridian  and  great 
circle  between  sub-vehicle  point  and  destination 
(Fig.  C-1).  [Symbol  3  also  used  for  longitude 
with  respect  to  vernal  equinox  line  (Fig.  1) J 

=  Flight  path  angle  (positive  above  horizontal) 
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Emissivity  of  vehicle  skin  at  stagnation  point. 

Heading  with  respect  to  north  (positive  clockwise). 
Vehicle  longitude  (positive  easterly) 

Longitude  of  tracking  station  i  (positive  easterly). 
Longitude  of  destination  (positive  easterly). 

Earth  central  angle  covered  by  re-entry  vehicle  at 
zero  bank  angle. 

Great  circle  arc  from  sub-vehicle  point  to  destination. 
(Fig.  C-1). 

Great  circle  arc  from  destination  to  reference 
great  circle.  (Fig.  C-1). 

Great  circle  arc  from  lateral  cardinal  point  of 
GAA  to  reference  great  circle.  (Fig.  C-1). 

Angle  between  vehicle* s  local  meridian  and  reference 
great  circle.  (Fig.  C-1). 

Local  air  density. 

Atmospheric  density  at  surface  of  the  earth. 

Bank  angle  (positive  for  right  wing  down). 

Vehicle  latitude  (positive  northerly) 

Latitude  of  tracking  station  i  (positive  northerly). 
Latitude  of  des ti nation  (positive  northerly) 

Angle  between  vehicle*s  initial  local  meridian 
and  great  circle  between  initial  sub-vehicle  point 
and  present  sub-vehicle  point. 

Angular  velocity  of  rotating  coordinate  system. 

Earth’s  angular  velocity. 
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